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Does the Universe Accelerate Equally in all Directions? 
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ABSTRACT 



We employ the Union compilation of Type la supernovae (SNe) with a maximum 
likelihood analysis to search for a dark energy dipole. To approach this problem, we 
present a simple, computationally efficient, and largely model independent method. 
We opted to weight each SN by its observed error estimate, so poorly measured SNe 
that deviate substantially from the Hubble law do not produce spurious results. We 
find, with very low significance, a dipole in the cosmic acceleration directed roughly 
towards the cosmic microwave background (CMB) dipole, but this is almost certainly 
coincidental. 
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1 INTRODUCTION 

The foundations for the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) cosmological model rest on the assump- 
tions of isotropy and homogeneity, both of which seem 
plausible on aver age. However, given the discovery of cos - 
mic acceleration dRiess et aZ.lll998l ; IPerlmutter et qZ.l ll.999). 
the FLRW model must attribute this observation to a cos- 
mological constant (dark energy) in Einstein's equations, 
which still lacks an acceptable p hysical interpretation (see 
iFrieman. Turner fe Hutererll2008l for a review). 

To alleviate these concerns, numerous authors have in- 
vestigated possible inhomogeneous, and void models t o ac- 
count for the o bserved acceleration (see Enqvist 2008| for a 
review, but also lAlnes. Amarzguioui fc G r0n 200(|and refer- 
ences therein). The typical approach is to employ the exact 
Lemaitre-Tolman-Bondi solutions for a radially inhomoge- 
neous, spherically symmetric Universe, and require that we 
must be almost centrally located in a large scale void. How- 
ever, a need to violate the Copernican principle has not yet 
been evid enced. Moreover, many of these models seem un- 
likely (see llshibashi fc Waldll2006l ). 

A more recent development by IWiltshirei (|2007al h sug- 
gests that cosmic acceleration can be explained in gen- 
eral relativity, by considering the differences in the quasi- 
local gravitational energy between observers in bound sys- 
tems (clusters of galaxies in bubble walls and filaments) 
and those in a freely expanding space (volume-averaged ob- 
servers), resulting in a cumulative time dilation effect due 
to the differences in their clock rates (see also IWiltshirei 
2007blU Using the Gold SN e sample l|Riess et aZ.I l200it ). 



provides a possible alternative to dark energy, although 
this model fails to provide a consistent fit to the Union 



Leith. Ng fc Wiltshire! (|2008l ) demonstrated that this model 



|Kowalski et aZ.ll2008r) and Constitution dHicken et alhoO^ ) 
SNe samples IjKwan. Francis fc Lewisll2009r ). contrary to the 
current concordant cosmology (ACDM), whereby all three 
SNe samples are fit by a consistent set of cosmological pa- 
rameters. 



Other authors have investig ated the possibility that 
dark energy is anisotropic ( see Coorav. Holz fc Caldwelll 
120081 and references therein). ICoorav et al\ l|2008h present 
a viable method to model dark energy inhomogeneities as 
a power spectrum of luminosity fluctuations. However, with 
such optimistic objectives, these studies will require a ded- 
icated all-sky survey to search for SNe. One can, how- 
ever, apply a simple method to measure the anisotropy of 
dark energy, or alternatively, search for a systematic direc- 
tional dependence in the SN e data l|Kolatt fc Lahavll200ll ; 
iGupta. Saini fc Laskarl [2008). This involves dividing the 
data into two hemispheres, and searching for the strongest 
directional dependence, which these authors find to have a 
significance of about la for previous compilations of SNe. 



* E-mail: rcooke@ast.cam.ac.uk; dlb@ast.cam.ac.uk 



In this article we employ the current SNe dataset and 
apply a new technique to test for the presence of a dark 
energy dipole. In §2 we detail the SN sample used in our 
analysis, and then formulate our dipolar model in §3. In §4 
we present the results of our dipolar modelling, and sum- 
marise our conclusions in §5. Throu ghout we ad o pt th e 
cosmological parameters derived by iDunklev et all |2009); 
H = 71.9 kms -1 Mpc" 1 , Q m = 0.258, and Q A = 0.742. 
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Figure 1. A projection of the spatial distribution of the Union 
SNe sample in Galactic coordinates. Note the relative uniformity 
of the points, except around the Galactic plane. The symbols 
correspond to those in Fig. [2] and are explained in £13.11 



2 THE TYPE IA SUPERNOVA SAMPLE 

For the present a nalysis, we employ th e Union sample of 
SNe, compiled bv lKowalski et all (|2008l fl. The Union com- 
pilation contains 307 identically processed SNe, drawn from 
a heterogeneous sample of 414. For th e present work, it is 
noteworthy that Kowal ski et al\ (|2008T ) corrected their data 
for the CMB dipole, and rejected SNe with CMB-centric 
redshiftifl z < 0.015. The Galactic coordinates for all SNe 
in the Union sample were obtained from the NASA/IPAC 
Extragalactic Database (NEDfl Indeed, as illustrated in 
Fig. [1] the Union SNe compilation are drawn from a spa- 
tially uniform sample, except of course the trivial lack of 
SNe observed in the Galactic plane. 

We acknowledg e the presence of a more recent sample 
of low redshift SNe ijHicken et alh oOV). and decide against 
including these additional data for the following two reasons. 
First, the signature of an anisotropic acceleration of the Uni- 
verse is only detectable at higher redshifts, where the effect 
of a cosmological constant is more exposed, and the motion 
of the local group is less of a contaminant. Therefore a low 
redshift sample should not affect the final result. Moreover, 
for this reason, we opted to include only the Union SNe 
with z > 0.2, leaving us with a total sample of 250. Second, 
the more recent data are processed using different routines, 
which could introduce a bias if we were to blindly combine 
both samples. With the goal of the present work borne in 
mind, we restrict our analysis to the Union compilation. 



3 ANALYSIS 

As outlined bv lKolatt fc Lahavl i|200ll ). a number of sources 
could produce deviations from isotropy in the SN data. From 
the results of their study, these authors raised concerns 
about the relatively small number of SNe, spread inhomoge- 
neously across the sky, with differing and uncertain errors. 
We therefore seek a method that is not biased by the clus- 
tering of the SNe data, whilst taking into consideration the 
individual uncertainties of each SN. Bearing in mind these 
characteristics, we employ a maximum likelihood strategy 



1 Available from: http : //supernova. lbl . gov/Union/ 

2 iKowalski et al\ j200Sf) find that their results do not depend 
significantly on the adopted redshift cutoff. 

3 http : //nedwww. ipac . caltech.edu/ 
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Figure 2. (a) Histogram of deviations from the Hubble law. (b) 
A Hubble plot of the entire Union SNe sample (grey symbols), 
and those with the largest deviations (crosses/circles coloured 
red/blue represent positive/negative deviation). The solid line 
corresponds to the Hubble law for the concordance cosmology. 



(unbiased to clustering), and devise a method that appro- 
priately weights the contribution of each SN based on its 
corresponding error measurement. 



3.1 The Hubble Deviation 

Before proceeding with a quantitative analysis of the Union 
sample, it is instructive to inspect whether any obvious pre- 
ferred deviations f rom the Hubble law exi s t in th e employed 
dataset. Following I Jha. Riess fc Kirshnerl (|2007T I, the devia- 
tion from the Hubble law is given by 



5H_ 
IT 



H d L - H d S f 



Ho dsN 

where dL is the luminosity distance in a k = Universe 
(l + z)c 



d L = 



Ho 



[n A + n M {i + z'yr 1/2 dz' 



(i) 



(2) 



and dsN can be calculated from the observed distance mod- 
ulus, [i — 51ogdsjv + 25. Fig. [5](a) illustrates that the data 
are only very slightly skewed; more SNe are above the Hub- 
ble law (i.e. with negative Hubble deviation) than below 
(positive Hubble deviation). Moreover, perhaps unsurpris- 
ingly, there are also a greater number of SNe that lie above 
the Hubble law when considering the subsample of SNe that 
exhibit the largest Hubble deviations ( \SH/H\ > 0.15), rep- 
resented in Fig. [1] and [2] by the crosses and open circles. 
The entire SNe sample are presented on a Hubble plot in 
Fig. 0(b). The fact that all of the largest deviations are 
exhibited by SNe at z > 0.2, validates our decision not to 
include SNe with z < 0.2. 
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The spatial distribution of the entire Union sample of 
SNe is presented in Fig. [T] where the symbols correspond to 
those in Fig. [2] Indeed, from this simple demonstration, an 
obvious dipole is not apparent. 



3.2 Maximum Likelihood Strategy 

We consider a k — Universe, and compare two models 
with and without an acceleration term. In essence, the result 
should not be affected by the model chosen to include the ac- 
celeration term, provided that the model represents the cen- 
troid of the data at a given redshift (i.e. selecting the ACDM 
model is sufficient). We define 5(A) = d(A) — d(A = 0) to 
be the luminosity distance deviation of an isotropically ac- 
celerating Universe, d(A), from a non-accelerating Universe, 
d(A — 0). The luminosity distance we would measure in a 
Universe with an anisotropic acceleration g!(Aa), can then 
be simply described by 



d(A A ) = d(A) + 5{A)X 



(3) 



where X varies over a range, say ±L, and provides a measure 
for the degree of anisotropic acceleration. For example, if 
X — 0, we simply recover d(A), and if X — —1, we are left 
with d(A = 0). Introducing a unit vector U that points in 
the direction of the i th SN, and a vector L, that points in the 
direction of 'preferred' acceleration, we can suitably choose 
X — U ■ L, thus Eq. 3 becomes 



5 M ,i = D(zi)(l + k ■ L) 
with, 

. _ d(A A ) - d(A = 0) 



D( Zi ) = 



d(A = 0) 

d(A) - d(A = 0) 
d(A = 0) 



(4) 

(5) 
(6) 



where D(zi) is the 'deviation parameter' for the i th SN with 
redshift z;. Using this formalism, we can introduce the mea- 
sured quantity Si (recall that d$N is the measured luminosity 
distance to a given SN), 



St 



d S N - d(A = 0) 
d(A = 0) 



(7) 



Thus, the vector L provides an estimate for the mag- 
nitude and direction of the dipole, whilst the deviation pa- 
rameter describes how the dipole changes with redshift. If 
a dipole in dark energy were to exist, we would expect it 
to be more exposed at higher redshifts where the accelera- 
tion due to dark energy is more apparent. We have therefore 
opted for a deviation parameter, and hence dipole, that in- 
creases with redshift. Ideally, one would prefer to divide the 
data into 'redshift shells' to derive an arbitrary form for the 
dipole, however, we are at present limited by the relatively 
small number of SNe. By choosing the model in this way, 
we are testing for the significance of a dipole in the accel- 
eration term in the SNe data. A dipole with a different z 
dependence would still be detected by our analysis but the 
significance level would not then be correctly estimated. 



3.2.1 The Likelihood Function 

The dipolar model (Eq. [4} forms the basis for the present 
analysis, and is employed using a maximum likelihood strat- 
egy. As the observations are not uniformly distributed 
around the sky, we pose the question, given the directions li 
of the observed SNe, what is the most probable vector L that 
gives an anisotropy of the form in Eq. [4p Non-uniform sky 
coverage will therefore only result in a larger error ellipsoid 
for L. 

Within this framework, poorly measured SNe that de- 
viate substantially from the concordance model will produce 
spurious results. We therefore weight each SN according to 
its observational error. Suppose the measurements Si come 
from true values Ai scattered by observational errors <7;. 
Then, given an observation Si, the probability of it arising 
from a true value Ai is 



Pr(Si\Ai) = 



2tto, 



exp 



-(Sj - A t f 
2a? 



(8) 



where Oi is the observers error estimate for the i th SN. Sim- 
ilarly, the probability that the true Ai arise from a scatter 
a about the proposed dipolar model (Eq. [4| , is 



Pr(A,\S A 



2na 



• exp 



(A, 



2a 2 



(9) 



The probability, therefore, of measuring the Si given our 
model is found by integrating over all true values 



(10) 



Pr(Si\S M ,i) = I Pr(5i|A0 Pr(A t \8 M ^ dA 4 

I / — (Si — SM,i) 2 



Pr(8i\5 M ,i) 



^2n(*i+af) eXP K 2(<r»+a?) 
The log-likelihood function is then given by 



£ = ln 



(11) 



3.2.2 Maximum Likelihood 

Substitution of the dipolar model into Eq. 1111 and maximis- 
ing the likelihood over all choices of the model parameters 
L and a, dC/dL = yields an equation of the form 

A-L = V (12) 

where A (a rank 2 tensor), and V are given by 

— D(z t ) 2 



A = 



■ li li 



v 



(Si - D( Zi ))D(zi) f 



dC/da — gives 



\ ^ 



SM,i) 2 



(13) 



(14) 



(15) 



We must now simultaneously solve the coupled non- 
linear Eqs. I12l and ll5l The most efficient pathway to a solu- 
tion is to construct a set of o values. Using Eq. 1121 one can 
invert A, and find L for each o. We then use Eq.[l5]to solve 
for the a value that maximises C. 
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Table 1. Parameter Fitting Results 



Model 


\L\ x 10 2 


«t 


bt 


cr x 10 2 


\L M \x 10 2 


£* 


L err /\L\ 


Pr(L) 


a Weighted Dipolar Model (z > 0.2) 


7.14 


287° 


46° 


0.00 


15.7 


105.5 


1.04 


18% (0.23cr) 


b Weighted Dipolar Model (0.2 < z < 0.56) 


3.90 


258° 


9° 


0.00 


19.1 


82.9 


3.44 


0.6 % (O.OO80-) 


c Weighted Dipolar Model (z > 0.56) 


13.8 


309° 


43° 


3.52 


26.1 


23.3 


0.83 


31 % (0.39<t) 


d Higher Quality List 


10.0 


118° 


25° 


14.2 


20.8 


97.1 


1.03 


18% (0.23<r) 



a,b,c,d Th ese models are referred to in the text as Sample 1,2,3,4. t Galactic Coordinates. ^Given by Eq. 1111 



3.2.3 Other Solution Pathways 

The non-linearity of the above equations is introduced 
through the weighting of each SN by its associated error. 
If we were to weight all SNe equally (i.e. en = constant for 
all i), Eqs. U3lT51 then take the form 



A = D(zi) 2 U k 

i 

V = J2( 5 i~ D{zi))D{zi)ii 

i 



(16) 
(17) 
(18) 



where N is the total number of SNe in the sample. We can 
now invert A and uniquely solve for L, and then use Eg. 1181 
to solve for the a that maximises the likelihood, C. 



4 RESULTS AND DISCUSSION 

In this section we present the results from our maximum 
likelihood analysis, and introduce a method to quantify the 
significance of our derived parameter values. Before proceed- 
ing, it is noteworthy that the adopted maximum likelihood 
strategy detailed in £13.21 requires all SNe coordinates to be 
converted from Galactic to Cartesian coordinate^. 



4.1 Significance of the Dipolar Models 

We ask ourselves, what is the probability that the vector 
L, that maximises the likelihood, is true? The covariance 
matrix for L is given by the inverse of Eq. [13] (i.e. A^ 1 ), or, 
for the solution pathway outlined in £13.2.31 a 2 AT 1 . Given 
the eigenvalues (Ai,A2,A3) and corresponding eigenvectors 
(ei, e.2, €3) of the covariance matrix, the three semi-principal 
axes of the error ellipsoid for L are si,2,3 = y Ai,2,3 ei,2,3- 
The density of probability associated with the vector L is 
therefore of the form 



p(L) = 



^(2^)3 X 1 \ 2 A 3 



,2 2 2 
/ e i , e 2 . e 3 

eXp -{2X- 1 + 2X- 2 + 2X- 3 



(19) 



where ei, e 2 and ez are the Cartesian coordinate axes in 
the direction of the principal axes of the error ellipsoid. The 
length of the vector L enclosed by the error ellipsoid is given 
by 

L err = [(L ■ si) 2 + (L ■ s 2 ) 2 + (L ■ s 3 ) 2 ] 1/2 (20) 



4 The x and z Cartesian axes are respectively directed towards 
the Galactic centre and North pole. 




Figure 3. An illustration of L and its error ellipsoid. The largest 
vector contained within L's error ellipsoid, Ljf, is also shown. 



We can now introduce the quantity fj, — L err /\L\, which de- 
fines the fraction of L enclosed by the error ellipsoid. Sup- 
pose we now expand (or contract) the error ellipsoid such 
that it contains a volume Ve, that just encompasses L. The 
probability that our derived L is true, is therefore given by 



Pr{L) = I p{L)dV E 

which after converting to spherical coordinates yields 



(21) 



Pr(L) = 



4ty 

(tt) 3 / 2 



u 2 exp y— u 2 ) du 



Pr{L) = eri(l/uV2) - J — exp ( - l/2^ 2 

" jJ, TV y 



(22) 



Additionally, one can find a vector, Lm, that is the 
largest vector contained within the error ellipsoid of L (see 
Fig. [3]). If we are certain that a dipole exists, and that L is 
within its error ellipsoid, the quantity \Lm\ is then a mea- 
sure of the maximum allowed increase in the acceleration, 
irrespective of the direction of L. 



4.2 The Weighted Dipolar Model 

We now consider the weighted dipolar model; weighting all 
SNe according to their measured errors, and implement the 
solution strategy outlined in § 13.2.21 The model parameters 
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Figure 4. For all figures, the y-axis corresponds to the case 
where dark energy is isotropic (i.e. X = li ■ L = 0), and the x-axis 
corresponds to the cosine of the angle between the i th SN and 
the 'preferred' direction L, weighted by the deviation parameter 
and the individual errors, so that poorly measured SNe contribute 
less, as required by the dipolar model. Thus, the dotted line has 
a gradient of \L\. For a perfect fit to the dipolar model, the data 
should lie along this dotted line. The data are averaged in 10 
equally spaced bins, where the errors represent the standard er- 
ror in the mean of the given bin. The domain is set to be the 
maximum abscissa of the unbinned data, (a) Weighted dipolar 
model (Sample 1) z > 0.2. (b) Weighted dipolar model (Sam- 
ple 2) 0.2 < z < 0.56. (c) Weighted dipolar model (Sample 3) 
z > 0.56. (d) Higher Quality List (Sample 4). 



Figure 5. For all figures, the y-axis again corresponds to the 
case where dark energy is isotropic (i.e. the same as Fig.|4jl, how- 
ever the x-axis now simply corresponds to the cosine of the angle 
between the i th SN and the 'preferred' direction L. Thus, we have 
used the derived L for each Sample, to determine whether or not 
there are any departures from isotropy in the weighted data. A 
general trend from bottom-left to top-right would be expected if a 
preferred direction did exist. The data are averaged in 10 equally 
spaced bins, where the errors represent the standard error in the 
mean of the given bin. (a) Weighted dipolar model (Sample 1) 
z > 0.2. (b) Weighted dipolar model (Sample 2) 0.2 < z < 0.56. 
(c) Weighted dipolar model (Sample 3) z > 0.56. (d) Higher Qual- 
ity List (Sample 4). 
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that maximise the likelihood function are presented in Ta- 
ble 1 (Sample 1), where the components of L are given in 
Galactic coordinates. For this case, a 2 is negligible compared 
to the minimum of in the sample (refer to Eq. 1 10 p . The fit- 
ting results are also illustrated in Fig. |4ja) , where the data 
are averaged in 10 equally spaced bins, and the error bars 
correspond to the standard error in the mean of each bin. 
Bins containing ^ 2 points are not shown. 

A striking observation of the results from the weighted 
dipolar model, is the similarity in the direction of L to 
the CMB dipole: (l,b ) = (263.99° ± 0.14°, 48.26° ±0.03°) 
(|Hinshaw et all 2009). Although there is a relatively mi- 
nor significance, this coincidence is certainly worth further 
scrutiny. 

To investigate whether the higher or lower redshift SNe, 
or both, favour this direction, we divide the SNe into two 
subsamples, each containing 125 SNe, and then run the 
weighted dipolar model routine on the two subsamples. This 
division corresponds to a redshift z — 0.56. The results are 
again presented in Table 1, and illustrated in Fig. Sfb) and 
2fc) for the two subsamples 0.2 < z < 0.56 (Sample 2) and 
2 > 0.56 (Sample 3) respectively. As can be appreciated 
from Fig. |3£b), and the final column of Table 1, the lower 
redshift subsample does not resemble the anisotropic model 
at all. The higher redshift subsample presented in Fig. |4jc), 
on the other hand, seems more consistent with the results 
of Sample 1. 

In order to test whether the Union compilation exhibit 
a dipolar acceleration, we plot the weighted data as a func- 
tion of Z; • L (i.e. the cosine of the angle between the i th SN 
and the 'preferred direction') in Fig. [5j a )~( c )- Should a pre- 
ferred direction exist, one would expect the data to show a 
general trend from the bottom-left quadrant to the top-right 
quadrant, a feature that is not observed in these samples. 
One could tentatively claim that this is true for Sample 3 
(in Fig. [5jc) ) , but the significance is very low. 



4.3 The Higher Quality List 

In § 13.2.31 an alternative linear method was described to 
solve for the model parameters, however, as mentioned pre- 
viously, one must be cautious when assigning each SN mea- 
surement an equal weighting, as poorly measured SNe that 
deviate substantially from the Hubble law could produce 
spurious results. These concerns should be alleviated by 
compiling a 'Higher Quality List', containing only the best 
quality SNe (Sample 4). We therefore restrict the sample 
to contain only the best measured SNe; SNe with uncer- 
tainties in magnitude > 0.35 are rejected. The distribution 
of accepted and rejected SNe are presented in Fig. H^a). 
The Hubble diagram for the Higher Quality List is shown in 

Fig.ETb). 

One would expect the results of the weighted dipolar 
model (Sample 1) to be consistent with the Higher Quality 
List, however, this does not seem to be the case. The direc- 
tion of L has changed by 108° , and therefore does not lend 
support to the results of Sample 1. For completeness, we also 
illustrate these results in Fig. UJd), and Fig. [SJd), however, 
neither demonstrate the trend that would be expected if an 
anisotropy existed in the Union SNe compilation. 
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z 

Figure 6. (a) Distribution of the SNe errors from the Union 
compilation. The light grey shaded regions correspond to the SNe 
with z < 0.2. The hatched regions correspond to the SNe that 
have Afi ^ 0.35. The unshaded histogram represents the Higher 
Quality List, with z > 0.2 and A/x < 0.35. (b) A Hubble plot 
of the Higher Quality List (symbols) and Hubble law for d(A) 
(upper solid line) and d(A = 0) (lower solid line). 

5 CONCLUSIONS 

The acceleration of the Universe appears to be equal in all 
directions. The anisotropy we find is a 14% ± 12 % increase 
toward (l,b) = (309°, 43°) and a corresponding decrease in 
the opposite direction (l,b) = (129°,— 43°), but the effect 
is only apparent in the higher redshift group with z > 0.56. 
These are those that exhibit the largest net acceleration. The 
direction of this weak anisotropy is only 31° from that of the 
CMB dipole as seen from the Sun, (l,b) = (264° , 48.3°), or 
17° from that of the dipole observed in the CMB frame 
of the elliptical g alaxies sample with v < 2000kms _1 , 
(l,b) = (311°, 26°) (|Lvnden-Bell et aZ.lll988D . The Union SN 
redshifts are corrected to the CMB system of rest, so their 

zero point cannot cause this eff ect. 

It is noteworth y that both lKolatt fc Lahavl l|200ll ) and 
iGupta et all (|2008l ) find a directional dependent systematic 
in the SNe data of similar significance to what we have found 
here, although, the directions are somewhat different. The 
present dataset cannot exclude a 26 % (i.e. 14+12) increase 
in the acceleration toward {l,b) — (309°, 43°) and together 
with a 26 % decrease in the opposite direction. More high 
accuracy z > 0.5 data from directions within 45° of (I, b) = 
(309°, 43°) or (l,b) = (129°, -43°) are needed to eliminate 
any possible anisotropy. 
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